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Turbulent channel flow of drag-reducing polymer solutions is simulated in minimal flow geometries. 
Even in the Newtonian limit, we find intervals of "hibernating" turbulence that display many features 
of the universal maximum drag reduction (MDR) asymptote observed in polymer solutions: weak 
streamwise vortices, nearly nonexistent streamwise variations and a mean velocity gradient that 
quantitatively matches experiments. As viscoelasticity increases, the frequency of these intervals 
also increases, while the intervals themselves are unchanged, leading to flows that increasingly 
resemble MDR. 



The energy dissipated in turbulent channel or pipe flow 
of a liquid can be dramatically reduced by low levels of 
long-chain polymer additives [IH3]. The most striking 
qualitative feature of this phenomenon is the existence 
of a so-called maximum drag reduction (MDR) asymp- 
tote pQ. For a given flow geometry at a given pres- 
sure drop (i.e. at a given Reynolds number Re), there is 
an asymptotic maximum flow rate that can be achieved 
through addition of polymers. Changing the concentra- 
tion, molecular weight or even the chemical structure of 
the additives has no effect on this asymptotic value. This 
universality is the major puzzle of drag reduction. 

Turbulent flow in or near the MDR regime displays 
important differences from Newtonian turbulence. Its 
most commonly discussed signature is a distinctive mean 
velocity profile U medin (y) that displays clear log-law be- 
havior well- approximated by a formula given by Virk: 
[/+ an = 11.7 In y+ - 17.0 pQ (Superscript "+" de- 
notes quantities nondimensionalized in inner velocity and 
length scales ^r w /p and r)/y/pr w ; r w is the time- and 
area averaged wall shear stress, r] and p are fluid viscosity 
and density and y is distance from the wall.) Addition- 
ally, streamwise vortices, which dominate near- wall dy- 
namics in Newtonian turbulence, are significantly weak- 
ened at MDR. Low-speed streaks become much less wavy 
in the streamwise direction and the streak spacing is sub- 
stantially larger @HZI- In addition, the Reynolds shear 
stress near MDR is substantially smaller than in Newto- 
nian turbulence [5HTT]. Indeed, the smallness of Reynolds 
shear stress is a central issue in a recent phenomenologi- 
cal model of MDR [12]. 

Based on these observations, many researchers have 
suggested that turbulence in the MDR regime is "transi- 
tional" [3] or "marginal" [12] in some sense that is not yet 
well-defined, and that the spatiotemporal flow structures 
that sustain turbulence in this regime are substantially 
different from those of normal Newtonian turbulence. In 
the latter case, one simulation approach that has been 
fruitful in identifying the self-sustaining structures is the 
so-called minimal flow unit (MFU) approach p~3j fT3] . 
This approach identifies the smallest flow domain (at a 



given Reynolds number) in which turbulence can be sus- 
tained. Accordingly, temporally intermittent phenomena 
can be identified more readily than in a large box, where 
spatial averages incorporate different regions in which 
the instantaneous behavior may be very different. The 
present work is the first to address the relevance of this 
approach for understanding MDR. 

We focus on plane Poiseuille flow with fixed pressure 
drop. Streamwise, wall-normal and spanwise directions 
are denoted x, y, and z, respectively. The no-slip bound- 
ary condition applies at y = zbl and periodic boundary 
conditions apply in x and z\ the periods in these direc- 
tions are L x and L z . Lengths are scaled with half-channel 
height / and velocities with the Newtonian laminar cen- 
terline velocity U at the given pressure drop. Time t is 
scaled with l/U and pressure p with pU 2 . The governing 
equations are: 

^-^„ + ^(v.. t ),v.„ = o (1 ) 

Wi ( Dot „ . „ , T \ b 

^(rrVt 1 -^)- (3) 

Eq. ([I]) describes conservation of momentum and mass, 
and the polymer conformation and stress tensors ol and 
t p are described by the FENE-P constitutive equation 
for a bead-spring dumbbell model of a polymer in so- 
lution ((Eqs. |2| and Q) [15 . Reynolds number Re = 
pUl/r], Weissenberg number Wi = 2XU/1 (A is the poly- 
mer relaxation time), viscosity ratio j3 = r] s /r] (r] s is the 
solvent contribution to zero-shear viscosity), and b is the 
upper limit of the extension of polymer. All simulations 
reported here are at Re = 3600, f3 = 0.97, b = 5000. The 
numerical integration procedure is described in [5]. 

A rigorous search for MFUs would consider the param- 
eter dependence of both and L+, a task involving an 



2 



(a) Newtonian 



<d v /9 y> 
x h 



bulk 



4- 
0.35 

3- 



u 1 u 



1000 2000 3000 

(b) l/l// =29 



4000 



<dv /dy>> 

d h 



bulk 



0.35 




<dv/dy>. 23^M 



1000 



2000 

t 



3000 



4000 



FIG. 1: Mean shear rates at the walls ("b"-bottom, "t"- 
top) and bulk velocity Uhu\k as functions of time for typical 
segments of (a) Newtonian and (b) viscoelastic ( Wi = 29) 
simulations. Solid and dashed lines show (dv x /dy) = 2 and 
(dv x /dy} = 1.8, respectively. In laminar flow, Uhuik = 2/3. 



impractically large number of simulation runs. Here we 
fix = 360 (which is in the range of streamwise sizes of 
Newtonian MFU Q]2[T4]) and vary L+ only. Although 
both length scales depend on parameters, Lf is arguably 
the quantity of more interest: the dominant structures 
are the streamwise streaks and the streamwise vortices 
aligned alongside them; restricts the streak spacing 
and the size of the vortices. 

Fig. [T] shows time series of instantaneous bulk average 
velocity Uh u \k and area-averaged shear rate (dv x /dy) at 
the top and bottom walls for (a) Newtonian flow and (b) 
viscoelastic flow at Wi — 29, where DR% = 26 and L+ 
has increased from 140 to 250. In the Newtonian case, 
Fig. [TJa), one occasionally observes long-lasting periods 
when the shear rate at one or both walls is substantially 
lower than the average value of 2 - for example the time 
interval 3200 < t < 3400. By momentum conservation, 
the bulk velocity increases during these periods. A simi- 
lar observation was made in the Newtonian MFU study 
of Webber et al. [14, . These periods will be termed "hi- 
bernation" , in contrast to the "active" turbulence found 
outside them. As Wi increases, hibernation periods be- 
come increasingly frequent (Fig. [TJb)) - since the bulk 
velocity increases during these periods, they contribute 
substantially to drag reduction. 

To systematically identify hibernation events, two cri- 
teria are used: (1) area- averaged wall shear rate at one or 
both walls drops below a cutoff value (dv x /dy)\ cu t E = 
1.8; and (2) it stays there for longer than a certain 
amount of time At cuto ff = 50. Hibernating periods so 
identified are shown in the middle panels of Fig. [I] as 
rectangular signals, on the top or bottom of the plot ac- 
cording to the wall(s) on which the criterion is satisfied. 

Based on this identification scheme, Fig. [5|a) shows, 
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FIG. 2: (a) Time scales and fraction of time spent in hiber- 
nation and (b) level of drag reduction and spanwise box size, 
as functions of Wi. (At the relatively low Reynolds number 
considered here, the flow laminarizes for Wi > 31). 
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FIG. 3: A hibernation event. Thick black lines are results 
at Wi = 29 (200 < t < 600 in Fig. [TJb)). Colored lines are 
from Newtonian simulations started at the dots, using velocity 
fields from the Wi = 29 simulation. 



as functions of Wi, the mean duration of the hiberna- 
tion periods T#, mean duration of active periods Ta, and 
fraction of time spent in hibernation Fjj- Corresponding 
results for minimal spanwise box size and level of drag 
reduction are in Fig. |2jb). Notice that the average du- 
ration Tjj of a hibernating period is almost completely 
insensitive to Wi. In contrast, the average duration Ta 
of an active turbulence phase decreases substantially af- 
ter onset of drag reduction. Therefore, at high Wi, vis- 
coelasticity compresses the lifetime of active turbulence 
intervals, while having virtually no effect on hibernation. 

The insensitivity of Th to Wi suggests that flow during 
hibernation does not strongly stretch polymer molecules. 
Indeed, at Wi = 29 the peak value of (a yy ), which 
is closely associated with streamwise vortex suppres- 
sion [I2j [I6j [17] , drops from about 210 in active turbu- 
lence to about 5 during hibernation, a 40-fold reduction. 
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FIG. 4: (a): Instantaneous mean velocity profiles of snap- 
shots i-v (colors), and time-averaged profiles in the Newtonian 
and Wi = 29 cases. (The latter are plotted in terms of con- 
ventional wall units.) (b): Instantaneous mean velocity pro- 
files from time instants iii and iv plotted in conventional wall 
units. For comparison, a downward- shifted plot of the Virk 
log-law is also shown, (c): flow structures of typical snap- 
shots in hibernation (iii) and active turbulence (v). Green 
sheets are isosurfaces v x = 0.3; pleats correspond to low- 
speed streaks; red tubes are isosurfaces of streamwise- vortex 
intensity Q2D = 0.02, calculated by applying the Q-criterion 
of vortex indentification [18] in the yz plane [5] [17] . 



These results suggest that hibernation should be very 
similar in the Newtonian and viscoelastic cases. To test 
this possibility, velocity fields from time instants before 
and during a hibernation event at Wi = 29 were used 
as initial conditions for a Newtonian simulation, the tra- 
jectories of which were then compared with those from 
the original viscoelastic simulation. Fig. [3] illustrates the 
original viscoelastic trajectory (thick black line) as well as 
Newtonian trajectories (colors) started at various times. 
For the Newtonian run starting before any sign of hi- 
bernation is observed (t = 205), active turbulence is sus- 
tained. However, runs started from later times show that 
once the system begins to enter hibernation, removing 
the polymer stress does not cause turbulence to revert 
to an active state, although the depth and duration of 
hibernation are weakly dependent on the start time. 

To better understand these results, we examine more 
closely the hibernating period shown in Fig. [3] Several 
time instants are selected as marked: (i) is just before 
turbulence enters hibernation; (ii) is on the path toward 



hibernation; (iii) and (iv) are within hibernation; (v) is 
after turbulence becomes reactivated. Fig. [4^i shows in- 
stantaneous area-averaged velocity profiles in the bottom 
half of the channel at these instants, plotted in inner units 
based on the instantaneous wall shear stress at the bot- 
tom wall (denoted by the superscript * rather than + ). 
By collapsing the viscous sublayer behavior onto a single 
curve, this choice of scaling best exposes the nature of 
and differences between the various time instants shown. 
In active turbulence (i and v), the profiles fluctuate sub- 
stantially. Profiles for snapshots completely in hiberna- 
tion (iii and iv) are fundamentally different. In particu- 
lar, in the range 15 < y* < 40, both profiles show a clear 
log-law relationship with a slope very close to (within 
10% of) the Virk MDR asymptotic slope of 11.7 pQ. 
(The differences between the intercepts is smaller than 
the scatter of the available data.) The figure also shows 
the time-averaged mean velocity profiles for the New- 
tonian flow (which has a log-law region in good agree- 
ment with the classical result /7+ ean = 2.51n?/ + + 5.5) 
and the flow at Wi = 29, which is, as expected, inter- 
mediate between the Newtonian and MDR profiles. To 
illustrate that the choice of scaling does not affect the 
conclusion that velocity profiles in hibernation display 
a log-law slope close to the Virk results, Fig. (4)3 shows 
the curves from time instants iii and iv in conventional 
wall units. Although shifted downward from the Virk 
profile (because in the instantaneous wall shear stress is 
less than the time-averaged wall shear stress), the log-law 
slope remains within 20% of the Virk value. The New- 
tonian hibernation periods (not shown, for brevity) are 
very similar and the Virk slope is observed there as well. 

Fig. [4^ shows flow structures corresponding to time 
instants iii and v. Within active periods (v), turbulence 
shows the expected highly 3D structure of streamwise 
vortices and low-speed streaks [13] [I9j [20] . During hiber- 
nation (iii), streamwise vortices are significantly weaker; 
low-speed streaks are still observed, but are weak and 
only weakly dependent on x. (The low shear events ob- 
served in the Newtonian MFU study of Webber et al [14] 
also display weak streamwise dependence.) Weak stream- 
wise vorticity and three-dimensionality are distinct char- 
acteristics of the MDR regime [U [3j SI d [7] . The weak 
effect of viscoelasticity on hibernating turbulence may 
lie in its nearly streamwise-invariant kinematics. In the 
limiting case of a streamwise invariant steady flow, mate- 
rial lines cannot stretch exponentially [21 ; accordingly, 
polymer stretch in such a flow will not be substantial. Fi- 
nally, the Reynolds shear stress during hibernation drops 
to very low values relative to active turbulence, where it 
peaks at about 0.8; the peak value during hibernation is 
about 0.3. Again, this result is consistent with observa- 
tions in the MDR regime [5HTT]. 

The qualitative picture that emerges from these simu- 
lations is thus the following. Active turbulence generates 
substantial stretching of polymer molecules. The result- 
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ing stresses act to suppress this turbulence and drive the 
flow toward a very weakly turbulent hibernating regime. 
During hibernation the polymer molecules are no longer 
strongly stretched and relax toward equilibrium. Eventu- 
ally, hibernation ends as new turbulent fluctuations begin 
to grow, and the system transits back into active tur- 
bulence. The active turbulence again stretches polymer 
chains and the (stochastic) cycle repeats. 

In this picture, experimental observations in which the 
Virk MDR mean velocity profile is found correspond to 
a limiting situation - not achieved at the low Reynolds 
number and small boxes studied here - where the frac- 
tion of time and space occupied by active turbulence be- 
comes small enough that the hibernating regime domi- 
nates the statistics. Active turbulence cannot vanish en- 
tirely, because it is known that on average, the polymer 
molecules carry a substantial fraction of the mean shear 
stress [9l [10] , and since hibernating turbulence does not 
stretch polymers, some active turbulence must remain. 
These considerations lead to a new picture of turbulence 
in the MDR regime as a state in which hibernating tur- 
bulence is the norm, with active turbulence arising inter- 
mittently in space and time only to be suppressed by the 
polymer stretching that it induces. 

This study focused on MFU flows at low Reynolds 
number, allowing analysis of an extensive data set in a 
regime where flow structures are relatively simple. Re- 
markably, even this regime displays clear signatures of 
the features of MDR commonly associated with higher 
Reynolds numbers. (Though it should be noted that the 
MDR asymptote is experimentally observed at values of 
Re all the way down into the laminar-turbulent transi- 
tion regime pQ.) To carefully evaluate the generality of 
picture just presented, future work will require simula- 
tions at higher Re and large boxes in combination with 
pattern analysis tools that can identify spatiotemporally 
localized regions of active and hibernating turbulence. 

In addition, attention must be focused on the hiberna- 
tion phenomenon. Waleffe has identified a class of non- 
linear traveling wave solutions in the plane Couette and 
Poiseuille geometries - saddle points in phase space - that 
share many characteristics with hibernating turbulence, 
specifically weak streamwise vortices and weak stream- 
wise dependence [22 . At least one family of these solu- 
tions has vanishing streamwise dependence and Reynolds 
shear stress as Re — >• oo. Similarly, other researchers 
[23H25] have identified both simple and chaotic saddle 
trajectories, "edge states", that lie on the boundary be- 
tween laminar and turbulent dynamics. These observa- 
tions may allow us to make rigorous the qualitative no- 
tion described in the introduction that the MDR regime 
is transitional or marginal - it may be that hibernating 
turbulence is very near in phase space to the boundary 
between laminar and turbulent flow. Finally, a better 
understanding of hibernation may lead to strategies for 



active flow control to maintain turbulence in a state of 
hibernation and thus dramatically reduce skin friction. 
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